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ABSTRACT 

Context. X-ray observations have unveiled the existence of a family of radio-quiet Isolated Neutron Stars whose X-ray emission 
is purely thermal, hence dubbed X-ray Dim Isolated Neutron Stars (XDINSs). While optical observations have allowed to 
relate the thermal emission to the neutron star cooling and to build the neutron star surface thermal map, IR observations are 
critical to pinpoint a spectral turnover produced by a so far unseen magnetospheric component, or by the presence of a fallback 
disk. The detection of such a turnover can provide further evidence of a link between this class of isolated neutron stars and 
the magnetars, which show a distinctive spectral flattening in the IR. 

Aims. We present the deepest IR observations ever of five XDINSs, which we use to constrain a spectral turnover in the IR 
and the presence of a fallback disk. 

Methods. We have used archived VLT observations of these neutron stars performed with the ISAAC instrument in the 
//-band (1.65 /im) and the available fallback disk models. 

Results. For none of our targets it was possible to identify the IR counterpart down to limiting magnitudes H ~ 21.5 — 22.9. 
Although these limits are the deepest ever obtained for neutron stars of this class, they are not deep enough to rule out the 
existence and the nature of a possible spectral flattening in the IR. We also derive, by using disk models, the upper limits on 
the mass inflow rate in a fallback disk. We find the existence of a putative fallback disk consistent (although not confirmed) 
with our observations. 
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1. Introduction 

Observations performed at 7-ray and X-ray energies have 
unveiled the existence of peculiar classes of Isolated 
Neutron Stars (INSs) which stand apart from the family 
of more classical radio pulsars as they are, with only one 
confirmed exception ( [Camilo et al.|2006 l, radio-quiet. The 
Soft Gamma-ray Repeaters (SGRs) and the Anomalous 
X-ray Pulsars (AXPs) are believed to be magnetars, i.e. 
young neutron stars with long spin periods (4-11 s) and 
hyper-strong magnetic fields {B ^ 10^^ — 10^^ G, see 
Woods fc Thompson|2006" I . The magnetar model (Duncan 
& Thompson 1992 1 would explain several observational 
properties of SGRs and AXPs, including the period, the 
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spin down rate and their luminosity, in excess to what is 
expected from a standard pulsar spin down, see [Woods fc:| 
Thompson (2006) for a recent review. Other radio-quiet 



INSs are the dim X-ray sources detected by ROSAT, here- 



after X-ray Dim INSs or XDINSs, see Haberl (20071 for a 



review. 

The XDINSs are characterized by steady emission pat- 
terns in the soft X-rays regime, where the spectra can be 
described by black-bodies with 40 < KT < 100 eV and 
by high X-ray to optical flux ratio {log{fx/ fopt) ~ 4 — 6). 
The XDINSs soft X-ray emission is likely powered by the 
cooling of the neutron star and it is pulsed with periods 
in the 3-11 s range. 



According to coohng models (Page et al. 20061, the 



inferred surface temperatures imply that XDINSs ages are 
in the range of ~ 10^ — 10^ years which, in the case of RX 
J1856.6-3754 and RX J0720.4-3125, is compatible with 
their place of origin being in the Upper Scorpio and in 
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Target Date Exp. time seeing airmass sky 

dd-mm-yyyy (s) (arcsec) 



RX J0420. 0-5022 


14-11-2003 


1980 


0.8 


1.4 


clear 




11-01-2004 


1980 


0.6 


1.3 


clouds 




14-01-2004 


1980 


1.1 


1.7 


cyrrus 


RX J0720.4-3125 


05-12-2003 


1980 


0.8 


1.0 


photometric 




09-01-2004 


1980 


0.6 


1.1 


cyrrus 




12-01-2004 


1980 


0.7 


1.0 


clear 


RX J0806.4-4123 


26-11-2004 


3960 


0.8 


1.0 


clear 




22-12-2004 


1980 


0.6 


1.1 


clear 


RX J1856.5-3754 


23-05-2003 


2760 


0.6 


1.1 


cyrrus 


RX J2143.0+0654 


07-06-2003 


2580 


0.5 


1.2 


cyrrus 



Table 1. Summary of the ISAAC observations of our targets. Column two gives the observing dates(s), columns 
three to five give the total integration time per observation, the average seeing and airmass. Column six gives the sky 
conditions. 



the Trumpler 10 OB associations respectively, once their 
present position and proper motion are taken into account 
( Kaplan et al][2007l ). 

There is evidence which suggest that magnet ars and 
XDINSs may be linked at some level. They have similar 
X-ray periods, and for four XDINSs, RX J1308.6-h2127 



(Haberl et al. 2003), RX J1605.3-H3249 (van Kerkwijk 



et al. 20041, RX J0720.4-3125 (Haberl et al. 20041 and 



RX J2143. 0-1-0654 (Zane et al 



20051 the observed X- 



ray absorption features, due either to proton cyclotron 
resonance or to atomic transitions, yield magnetic fields 
B ~ 6 - 7 X 1013 - 1.0 X 10"G. For RX J0720.4-3125 
and RX J1308. 6-1-2127 the measured period derivatives 
■P ~ 7 X lO-i'i s s-i an d P - 11 x lO^^* s s'^ ( [Kaplan fc 
van Kerkwijk 2005a|b ^ imply, assuming magnetic dipole 
spin down, S ~ 3 x 10^'^ G, i.e. consistent with the values 
derived from X-ray spectroscopy. In addition, for the mea- 



sured parallactic distance of 360 pc (Kaplan et al. 20071 



the X-ray luminosity of RX J0720. 4-3125 turns out to be 
larger than its inferred spin down energy, a characteristic 
which makes it similar to the magnetars which have an X- 
ray luminosity by at least 2 orders of magnitudes higher 
with respect to their spin down energy. In the IR, magne- 
tars have a peculiar flux distribution which flattens with 
respect to the extrapolation of the X-ray spectra, see e.g. 



Israel et al. (20051. This flattening could be due either to 



the contribution of emission from a fallback disk, or to a 
turnover in the magnetospheric emission from the magne- 
tar. In turn, the emission from the fallback disk could be 
generated by mass transfer and viscous dissipation within 
the disk and/or by re-processing of the X-rays from the 
neutron star. Were the XDINSs show a similar energy dis- 
tribution as the magnetars, flattening their spectra in the 
IR, it would strengthen the link between these two classes 
of objects. 

In this paper we report on recent IR observations of 
five XDINSs: RX J0420.0-5022, RX J0720.4-3125, RX 
J0806.4-4123, RX J1856.6-3754 and RX J2143. 6-^0654. 
The paper is divided as follows: observations, data reduc- 
tion and results are described in §2, while the discussion 
is presented in §3. 



2. Observations, Data Reduction and Results 

XDINSs IR observations were conducted using the NIR 



spectra-imager ISAAC (Moorwood et al. 19981 mounted 
at the First Unit Telescope (UTl) of the ESQ Very 
Large Telescope (VLT) at the Paranal Observatory in 
Chile under program IDs 071.C-0189(A), 072.C-0051(A) 
and 074.C-0596(A), Pis R. Neuhauser and B. Posselt. 
The instrument configuration was set to the Short 
Wavelength (SW) camera, equipped with a Rockwell 
Hawaii 1024x 1024 pixel Hg:Cd:Te array, which has a pixel 
size of 0'.'148 and a field of view of 152x152 arcsec. All 
observations were performed through the 7J-band filter 
(A = 1.65/im; AA = 0.30/im). To allow for subtraction of 
the variable IR sky background, each exposure was split 
in sequences of 14 to 33 short randomly dithered expo- 
sures with Detector Integration Times (DIT) of 12 s, and 
5 NDIT repetitions along each point of the dithering pat- 
tern. Observations were performed between May 2003 and 
December 2004. Table [T] reports the observations sum- 
mary. Science exposures were retrieved through the public 
ESO archiveQto gether with the closest in time calibration 
files. 

T he data were re duced using the ESO's echpse pack- 
age (Devillard 2001 IP for de-jitter and sky subtraction. 



Unfortunately, we could not find in the ESO archive IR 
standard stars observed sufficiently close in date and time 
to our targets. For this reason, we have used the H- 
band fiux of a number of 2MASS stars identified in the 
frames as relative calibration sources. Since the photomet- 
ric accuracy of 2MASS is as good as ~ 0.1 magnitudes, 
and since a direct on-the-frame photometric calibration 
is not affected by overnight sky background fluctuations, 
we conclude that the accuracy of our photometric cali- 
bration is fully satisfactory. We have used all available 
2MASS stars in each frame (between 3 and 7) for the rel- 
ative photometric calibration. Also, the wavelength over- 
lap of the 2MASS and the ISAAC H-hand filters is almost 
complete and the introduced error is within 0.03 magni- 



^ http://archive.eso.org/ 
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Target 




R.A. 


Dec. 


A r 


rcf. 


RX J0420.0- 


-5022 


04 20 01.95 


-50 22 48.1 


0'.'6 


1 


RX J0720.4- 


-3125 


07 20 24.96 


-31 25 50.1 


0'.'2 


2 


RX J0806.4- 


-4123 


08 06 23.40 


-41 22 30.9 


0'.'6 


1 


RX J1856.5- 


-3754 


18 56 35.62 


-37 54 35.3 


0'.'2 


3 


RX J2143.0+0654 


21 43 03.38 


+ 06 54 17.5 


0'.'6 


4 



Rea et al. (20071. 



& Kaplan (2006), 3: van Kerkwijk & Kulkarni (20011, 4 



tudes. Photometry on the 2MASS stars was performed via 
standard aperture photometry with the Skycat - GAIA 



package (P.W. Draper & Berry 2001). Frame distortions 



(amounting up to 2 pixels and depending on position of 
the source in the frame) were not corrected for, and intro- 
duce an error of up to 0.1 magnitudes on the photometry 
of the individual 2MASS stars, which largely dominates 
over the statistical errors. We finally derive a global pho- 
tometric accuracy of ~ 0.1 — 0.2 magnitudes. 

The astrometric solution for the stacked ISAAC im- 
ages was computed by using as a reference the positions 
and coordinates of stars selected from the 2MASS cata- 
logue, yielding an rms between 0'.'09 and 0'.'14. To regis- 
ter the positions of our targets we have used as reference 
the most precise coordinates available in the literature 
(see Table ^. For RX J0420.0-5022 and RX J0806.4- 
4123 we have used the X-ray coordinates obtained with 



Chandra (Haberl et al. 20041, for RX J0720.4-3125 and 



RX J1856. 6-3754 we have used the coordinates of the op- 
tical counterparts after correcting for the measured proper 



motions ( 


van Kerkwijk & Kaplan 


12006 and 


van Kerkwijk 


& Kulkarni 2001 respectively) while for RX J2143.0-F0654 



Table 2. Coordinates and position uncertainties of the 
observed objects. The uncertainties refer to a 90% con- 
fidence level for the sources not detected in optical (RX 
J0420.0-5022, RX J0806.4-4123, RX J2143.0-h0654). The 
reference epochs for the positions are detailed in the ref- 
erenced papers: 1: Haberl et al. (20041, 2: van Kerkwijk 




10 arcseo. .>'■ 



Fig. 1. ISAAC iJ-band images of the field of RX 
J2143. 6-1-0654. Although two sources appear within the 
large XMM error circle, no source is visible within the 
smaller Chandra error circle. Moreover the two sources 
within the XMM circle have large values of B-V, not con- 
sistent with the nature of neutron stars emission, and have 
optical luminosity a factor 10 higher than expected (see 
Mignani et al.||2007al and the text). 



NS Name 
RX 



J0420.0-5022 
J0720.4-3125 
J0806.4-4123 
J1856.5-3754 
J2143.0-H0654 



H 

magnitude 



> 21.7 

> 22.7 

> 22.9 

> 21.5 

> 22.1 



d 
pc 



10 



« 345 
3601^^° 
« 250 

« 430 



Nh 
20 cm-2 


A„ 
mag. 


1.6 


0.016 


1.2 


0.012 


1.0 


0.010 


0.7 


0.007 


2.4 


0.024 



Table 3. Summary of the iJ-band flux upper limits (ex- 
pressed in magnitude) for the five XDINSs. Limiting mag- 
nitudes are computed for a 3(T detection level. Columns 
three and four give the neutron star distance and the Nh, 
column five gives the interstellar extinction in the TJ-band, 
Ah, computed from the Nh derived from the X-ray spec- 



we use the coordinates obtained with XMM (Zane et al. tral fits using the relation of Paresce (19841 as done by, e.g. 



2005). The overall uncertainties on the targets position 
account for the errors on the sky coordinates (see Table 

for the rms of the astrometric solution, for the abso- 
lute uncertainty of the 2MASS coordinates (^ 0'.'2) and, 
for RX J0720.4-3125 and RX J1856.6~3754, for the errors 
on the proper motion propagation between the reference 
and the observation epochs. 

We do not detect any candidate IR counterpart for four 
of our sources, even extending the search area to 3 error 
circle radii (confirmed independently in Posselt|[2007 1. In 
the case of RX J2143. 0-1-0654, four possible counterparts 



Posselt et al.| ( |2007[ ). We refer to [Posselt et al.| ( |2007[ ) for 
the estimation of Nh, and the distances of RX J0420.0- 
502, RX J0806.4-4123 and RX J2143.0-h0654. For the ob- 
jects with measured parallaxes, RX J0720. 4-3125 and RX 
J1856. 5-3754, the distance measurements are taken from 



Kaplan et al. (2007) and from van Kerkwijk & Kaplan 



(2006), respectively. 



were already identified in Mignani et al. (2007a) based 



on their proximity with the XMM large error circle (Fig. 
1), but they were all considered unlikely candidates due to 
their large B—V > 0.5, to their optical brightness, a factor 
10 higher than expected and, at least for two objects, to 
their apparent extension; the only point-like object within 
the XMM error circle in the IR frame has a low detection 



significance (< 5a). Very recently the Chandra position 
was published for this object ( Rea et al.|2007 ). No sources 
are visible within the smaller Chandra error circle. 

The derived upper limits, at the 3<t level, are shown in 
table m 
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3. Discussion 

We note that our derived IR flux upper limits are well 
above the extrapolation of the XDINSs' black body X- 



ray spectra, e.g. Haberl (20071, which predicts values of 



H > 25. For RX J0720-3125 and RX J1856-3754, the 
only two XDINSs for which a reliable characterization 
of the optical flux distribution has been obtained 



[Kerkwijk fc Kulkarni^ 2001} |Motch et al.||2003"| , an hypo- 
thetical spectral flattening redward of the i?-band would 
imply H < 24 and H < 22, respectively. The estimation 
is made assuming a constant flux normalized to the mea- 
sured i?-band magnitude of the targets. Comparable val- 
ues can be expected for the other XDINSs in our sample 
assuming similar optical spectra and and X-ray to optical 
flux normalizations. 

We have used the measured _ff-band flux upper limits 
of the five XDINSs to derive their corresponding IR lumi- 
nosities upper limits. Fluxes have been normalized either 



for the measured optical parallax distances ( van Kerkwijk 



& Kaplan||2006 


Kaplan et al. 


20071 


distances presented in Posselt et al. 



interstellar extinction correction has been computed ac- 



cording to the relation of Fitzpatrick ( 1999 ) using as a 
reference the Ay inferred from the derived from the 



X-ray spectral fits in Posselt et al. (20071 and the rela- 
tions of Predehl & Schmitt (1995) and |Paresce] ( |1984[ ). In 
both cases, the computed -ff-band interstellar extinction 
is significantly smaller than the global photometric uncer- 
tainty. 

We have compared the XDINSs' IR luminosity upper 
limits with the measured IR luminosities of both rotation- 
powered pulsars and magnetars. For these objects, IR lu- 
minosities have been computed from the measured and 



the extrapolated K-hand magnitudes (see Mignani et al 



2007b for details) using the same procedure outlined 
above. Although the comparison between IR luminosities 
obtained in different bands is not formally correct, the 
H and i^-bands are close enough in wavelengths so that 
the error due to the unknown color correction is negligible 
with respect to the overall luminosity error budget. The 
comparison between different INS classes is shown in Fig. 
2, where IR luminosities are plotted vs. the neutron stars' 
rotational energy loss E (assuming / = 10'*^ g cm^). As 
pointed out by Mignani et al. (2007b I, the correlation be- 



tween the IR luminosity of rotation-powered pulsars and 
E suggests that it is of magnetospheric origin, while the 
intrinsically higher IR luminosity of the magnetars is likely 
powered by the emission from a fallback disk and / or by the 
magnetic field decay. Since only upper limits on Ljji and 
E (with the only exception of RX J0720-3125) are avail- 
able for the XDINSs, it is impossible to recognize a pattern 
in the diagram. XDINSs could either follow the same Lju 
vs E correlation of rotation-powered pulsars (or alike), 
and thus be magnetospheric emitter in the IR, or they 
might not follow any correlation at all, and thus be, e.g. 
mostly thermal emitters, as they are in the optical and X- 
rays. Alternatively, they could represent an intermediate 




30 32 34 36 38 
Log(E) (erg/s) 



Fig. 2. IR luminosities of isolated NSs and derived 
XDINSs upper limits plo tted as a function of the NSs 

updated 



rotational energy loss E (Mignani et al 



2007b 



from). Empty circles, filled squares and filled triangles in- 
dicate rotation-powered pulsars, magnetars, and XDINSs 
respectively. The dashed line corresponds to the linear fit 
for rotation-powered pulsars while the solid line shows the 
limit case Log(Ljfi) = Log{E). For the rotation-powered 
pulsars and the magnetars the luminosity has been com- 
puted in the K-hand, while for the XDINSs the luminosity 
upper limits refer to the _ff-band. 



class between the magnetars and the old rotation-powered 
pulsars, where all emission processes may coexist. In any 
case, their low IR emission clearly segregates them from 
magnetars. Whether this is due to an intrinsic difference 
between the two classes, or it is simply due to the lower 
X-ray luminosity of the XDINSs which would reduce the 
contribution of the putative disk IR emission, or to a low 
mass inflow rate through such a disk (or both) is presently 
unclear. 

Searches for fallback disks in different families of neu- 
tron stars were recently undertaken by several authors, 



e.g. Wang et al. (2007), Mignani et al. (2007b I on four 



compact central objects and on one high magnetic field 
radio pulsar respectively, both reporting non-detections. 



Instead Wang et al. ( 2006 1 successfully detected a fiatten- 
ing in the spectrum of the AXP 4U 0142-1-61 from Keck 
(fC-band) and Spitzer observations, interpreting the data 
as a detection of a fallback disk of mass « 1O~^M0 around 
the neutron star. The colours of the detected source are 
unusually red with respect to main sequence stars or giant 
stars. Fallback disks could be found in principle around 



XDINSs as well (Perna et al. 2000 Michel & Dessler 1981 1. 



To set constraints on the presence of underlying fall- 
back disks around the XDINSs in our sample we have 
used the derived _ff-band flux upper limits. Since these 
objects do not display radio emission, an hypothetical 
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Target 


Lx 


P 


P 


M {Rin — Rcor) 


M {Rir, = Rlc) 


Ref. 




10^^ erg s"^ 


s 


10-"ss-i 


IO'^^Mq year"^ 


1O"^°M0 year"^ 




RX J0420. 0-5022 


3.2 


3.45 


< 92 


0.56 


2.06 


1 


RX J0720.4-3125 


33.7 


8.39 


0.698 




3.17 


2 


RX J0806.4-4123 


3.5 


11.37 


< 18 


0.48 


3.49 


1 


RX J1856.5-3754 


2.8 


7.06 


< 19 


0.40 


2.70 


3 


RX J2143.0+0654 


20.3 


9.44 


< 60 


0.67 


3.81 


4 



Table 4. Estimates of the upper limits of the accretion from a fallback disk. The luminosities reported in column two 
are referred to the distances in column three of table [3] In column 5 and 6 are shown the upper limits of the accretion 
rates (in units of \Q~^^Mq year~^) for the inner disk radius equal to the corotation radius and to the light cylinder 
radius, respectively. The X ray luminosities are obtained from the measured XMM fluxes in the 0.1-2.4 keV energy 



band ( 


Haberl||2004| 


Zane et al. 


( 


2005 


) 



disk could extend well inside the light cylinder radius. We 
note that the detection of pulsed radio emission at very 
low frequencies (42-112 MHz) has been recently reported 
for RX J13 08.6-I-2127 and RX J2143.0-h0654 (Malofeev^ 
et al.||2007 l but has not been independently confirmed so 
far. The inner radius of the disk is generally expected at 
the magnetospheric radius. Since this is dependent on the 
magnetic field and on the accretion rate of the disk which 
is not known a priori, this radius cannot be computed 



( Chatterjee et alj2000 ). However, given the very low X-ray 



luminosity level of the XDINSs, we do not expect that an 
hypothetical fallback disk would be accreting on the neu- 
tron star, which implies that the minimum value that the 
magnetospheric radius can take is the corotation radius. 
Since the corotation radius only depends on the neutron 
star period, it can be computed for all of our objects. The 
disk emission was computed using the models of |Perna| 
et al. (2000|, which include both the contribution to the 



emission from viscous dissipation, as well as the contribu- 
tion from reprocessing of the X-ray fiux. However, given 
the very low X-ray luminosity level of most XDINSs (see 
Table 4), we found that, for fiuxes below (but close) to the 
IR limits, the contribution to the disk emission is generally 
dominated by the viscous dissipation part. Therefore, our 
observational limits translate into limits on the maximum 
possible mass infiow rate through the disk. The exception 
is represented by RX J0720. 4-3125, which has the largest 
X-ray luminosity and for which the fiux from a disk with 
inner edge at the co-rotation radius would be instead dom- 
inated by the reprocessing of the X-ray emission from the 
neutron star. The expected disk IR emission would in this 
case be brighter than our derived upper limit, unless the 
disk inner radius is much larger than the co-rotation ra- 
dius or the disk has an inclination with respect to the line 
of sight higher than the assumed 60°. 

The results are shown in Table [4] where the maximum 
M compatible with the limits is reported. For each ob- 
ject, we derived the limits using two different values for 
the inner radius of the disk: = Rcor and R^ = Ric, 
where Rcor is the corotation radius (the absolute min- 
imum), and Rlc is the light cylinder radius (for refer- 
ence). Clearly, when the radius is closer in (i.e. when 



Rin — Rcor)i the limits on M are more stringent, due 
to the fact that there is a larger contribution of fiux com- 
ing from the inner regions of the disk. This translates in 
an upper limit of the mass infiow though the disk of at 
least « 0.4 • lO^^^Mgyear-i. 

4. Conclusions 

We report on the deepest IR observations ever of 5 of the 
7 known XDINSs. None of the sources was detected in 
our data at the 3a level. Limiting iJ-band magnitudes are 
in the « 21.5 — 22.9 range. These limits do not allow to 
either exclude or confirm a spectral fiattening red-ward of 
the i?-band as observed in some magnetars. 

We have compared the derived IR luminosity upper 
limits with the measured values of other classes of isolated 
neutron stars. Interestingly, we found that the distribution 
of our targets in the L//j — E plane shows a separation 
between the XDINSs and the magnetars due to their lower 
(by at least two orders of magnitude) IR luminosity. Due 
to the indetermination we have on the IR fiux and on P 
we cannot exclude nor confirm that XDINSs will follow a 
trend similar to the rotation-powered pulsars in the Lj — 
E plane, implying IR emission from the magnetosphere. 

We have investigated the possible existence of fallback 
disks around our targets and found that the contribution 
of the disk emission would be mostly dominated by viscous 
dissipation. Although we cannot confirm the existence of 
fallback disks, we have derived disk mass infiow rate upper 
limits which are consistent with the age and the disk mass 
estimates of our targets. 

Multi-color deeper IR observations, both in the near 
infrared and in the medium infrared (Spitzer) bands, could 
better probe the possibility of spectral fiattening in the IR 
and its origin. 
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